Phase shifts in nonresonant coherent excitation 
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Far-off-resonant pulsed laser fields produce negligible excitation between two atomic states but 
may induce considerable phase shifts. The acquired phases are usually calculated by using the 
adiabatic-elimination approximation. We analyze the accuracy of this approximation and derive 
the conditions for its applicability to the calculation of the phases. We account for various sources 
of imperfections, ranging from higher terms in the adiabatic-elimination expansion and irreversible 
population loss to couplings to additional states. We find that, as far as the phase shifts are 
concerned, the adiabatic elimination is accurate only for a very large detuning. We show that the 
adiabatic approximation is a far more accurate method for evaluating the phase shifts, with a vast 
domain of validity; the accuracy is further enhanced by superadiabatic corrections, which reduce the 
error well below 10 -4 . Moreover, owing to the effect of adiabatic population return, the adiabatic 
and superadiabatic approximations allow one to calculate the phase shifts even for a moderately 
large detuning, and even when the peak Rabi frequency is larger than the detuning; in these regimes 
the adiabatic elimination is completely inapplicable. We also derive several exact expressions for 
the phases using exactly soluble two-state and three-state analytical models. 

PACS numbers: 03.65. Vf, 03.67.Ac, 42.50.Dv, 03.67.Bg 



I. INTRODUCTION 

Far-off-resonant laser pulses are a popular tool for in- 
ducing controllable phase shifts in atomic states. These 
phase shifts — usually referred to as dynamic Stark shifts 
- are frequently used in the construction of dynamic 
phase gates, which are a basic tool in quantum informa- 
tion processing [1]. In many algorithms (e.g. Grover's 
quantum search [2] ) one has to prepare such phase shifts 
very accurately [3]. Insofar as quantum algorithms in- 
volve a great number of phase gates, the accuracy of the 
latter is of crucial importance for high-fidelity quantum 
information processing. 

There are three major types of phase gates: dynamic 
[3], geometric [4] and using relative laser phases [5]. 
While the latter two types have certain advantages in 
terms of robustness against parameter fluctuations, these 
come at the cost of more demanding implementations. 
The dynamic phase gate benefits from the simplicity of 
implementation (because, unlike the other phase gates, 
it requires just a single off-resonant pulsed field), which 
determines its wide-spread use. 

We emphasize that such phase shifts also emerge in 
various more traditional dynamical problems involving 
complicated linkage patterns. The quantum dynamics of 
these multistate systems can often be understood only by 
reduction to simpler two- or three-state systems by us- 
ing adiabatic elimination of all far-off-resonant (virtual) 
states. For instance, it is mandatory to account for such 
dynamic Stark shifts in excitation of multiphoton transi- 
tions by femtosecond laser pulses [6]. 

As far as a phase shift of ir is concerned the simplest 
approach is to use a resonant 2tt pulse. A variable phase 
shift (p, however, requires a field with a suitable detuning 
and intensity; such a variable phase shift is required, for 



example, in deterministic quantum search [7]. 

The dynamic Stark shift, and the ensuing phase shift, 
are usually calculated by eliminating adiabatically the 
off-resonant state(s). In this paper we show that, unless 
applied very carefully, the adiabatic elimination (AE) ap- 
proximation can lead to significant errors in the value of 
the phase. We analyze various sources of errors, ranging 
from higher terms in the AE expansion to population de- 
cay and shifts from additional states [8], and show that 
the standard AE approximation is accurate only for very 
large detuning A (Sec. III). We then present a method 
for the evaluation of the phase based on the adiabatic 
approximation (Sec. IV). This approximation provides 
a simple formula for the gate phase, which contains the 
AE phase as a limiting case for |A| — > oo, but it is also 
valid for moderately large detunings (|A|T > 1, with T 
being the pulse width). Then we include superadiabatic 
corrections (Sec. IV C), dissipation (Sec. IV D), effects 
of additional states (Sec. V), and several exact solutions 
(Sec. VI). 



II. BACKGROUND 

We consider a two-state quantum system (a qubit) in- 
teracting with a coherent field (Fig. I, left) and we wish 
to estimate the accumulated phase during this interac- 
tion. The phase gate is defined as 
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where F = U(£f) is the desired form of the propagator 
U(i) at the final time if. The propagator U(i) satisfies 
the Schrodinger equation [9], 



ihd t U(t) = H(f)U(i), 



(2) 
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FIG. 1: Linkage diagram for a two-state system (left), a three- 
state system in a ladder configuration (middle), and a three- 
state system in a V configuration (right). 



with the initial condition U(ti) = I at time t\. This prop- 
agator allows us to calculate the phases of the two states, 
accumulated during the interaction, for arbitrary initial 
conditions. The Hamiltonian in the interaction represen- 
tation and the rotating-wave approximation reads 



H(t) 



o n(t) e - <z, w 

n(t)e iD W 



(3) 



with D — A(i') dt', where A = ujq — uo is the detun- 
ing between the laser carrier frequency u> and the Bohr 
transition frequency ujq. For simplicity, we assume here- 
after that the detuning is constant, A = const; the re- 
sults can be readily extended to time-dependent A(t). 
The Rabi frequency Cl(t) — — d.E(i)/7i parameterizes the 
coupling between the electric field with an envelope E(t) 
and the transition dipole moment d of the system. 

In terms of the two probability amplitudes of the qubit 
states c\(t) and c 2 (t), the Schrodingcr equation reads 



ihd t c(t) = H(t)c(i), 



(4) 



with c(t) — [ci(t), C2(t)] T . For initial conditions Ci(U) = 
1 and c 2 (ti) = 0, the action of the phase gate F reads 



Fci(t i ) = e , *ci(t i ) ) Fc2(ti)=0. 



(5) 



For initial conditions ci(U) = and c 2 (t\) = 1, we have 
Fci(ii) = 0, Fc 2 (ii) = e-^c 2 (ti). (6) 

Solving the Schrodingcr equation (4) for the initial condi- 
tions Ci(ti) = 1, c 2 (ij) = is sufficient for the calculation 
of the gate phase: if the Hamiltonian (3) produces the 
phase change (5), it will produce the phase gate (1). 



III. ADIABATIC-ELIMINATION 
APPROXIMATION 

A. Steady-state solution 

We begin with the traditional adiabatic elimination. 
To this end, it is suitable, with the phase transforma- 
tion d(t) = h(t), c 2 (t) = b 2 {t)e- lD{ - t \ to write the 



Schrodinger equation (4) in the energy picture, 



idMt) 
id t b 2 (t) 



£fi(t)6i(t) + A&2(t) 



(7a) 
(7b) 



with the initial conditions b\(ti) = I and b 2 (t[) — 0. The 
AE approximation is applicable when the field is tuned 
far off resonance (|A| ^> O), which implies a small tran- 
sition probability. Then we set b 2 (t) = 0, find b 2 (t) from 
Eq. (7b), and substitute it in Eq. (7a). The result is a so- 
lution of the form (5), with a zero transition probability 
and a phase factor 



i 



O(t) 2 
4A 



dt. 



(8) 



In fact, this is the steady-state solution, that is just the 
first term in an asymptotic expansion over A. Unfortu- 
nately, this expression has only a small region of validity 
|A| » Q and does not give us a rigorous error estimation. 



B. Adiabatic elimination: higher terms 

In order to find the next terms in the asymptotic ex- 
pansion over A it is more convenient to start from the 
original interaction representation, Eqs. (3) and (4), 



id tCl (t) = \n{t)e- iD ^c 2 {t), 
id t c 2 (t) = \n(ty D{t ^ Cl {t). 



(9a) 
(9b) 



A formal integration gives 



C2(t) 



n{t'y At ' Cl {t')dt'. (fo) 



Now we integrate by parts, 



C2 



(t) - -/_ 



* Cll+I 



^gW)de^' 



Cl(t)e + / ^\ dv di' 

J — OO 



n(t) 

' 2A ' 

^ ClWeiAt + i 9t[fi(t)clW]e4At 



1 



2iA 2 



e* At 'd t , t , m') Cl (t')} dt'. 



We substitute this expression in Eq. (9a) and obtain 

Cl (i) = e -7(*) e ^W, (12) 

where 7 w for a large detuning. Indeed, for smooth 
pulse shapes, the transition probability (1 — e -27 here) 
vanishes exponentially with A. For example, for a 
Gaussian pulse, the transition probability vanishes as 
~ sech 2 [7rAT/21n(f7 T)] [10], whereas for a hyperbolic- 
secant pulse it vanishes as ~ sech 2 (7rAT/2) [11]. 
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The phase in Eq. (12) reads 



4A 



16A 3 



o 6 

A 5 



■Ol-rr) dt. 



The calculation of the higher terms is increasingly com- 
plicated and barely useful. We will refer to this phase as 
AE2, in order to distinguish it from the AE phase (8). 
This expression shows that the AE approximation is good 
only for a very large detuning and a smooth pulse (be- 
cause of the presence of Q(t)). The reason is that in an 
expansion for a phase it is not sufficient to retain the lead- 
ing term and demand the next term to be much smaller 
than it; one must also demand all neglected terms to be 
much smaller than unity. Obviously, terms of the order 
of unity or larger cannot be discarded in the estimate of 
the phase because the latter is defined modulo 2tt; hence 
the leading term alone may provide a value that is not 
even close to the exact value. We emphasize, however, 
that if the transition probability is concerned, then the 
leading term in the AE approximation, which is of order 
0(fll/A 2 ), provides an adequate estimate. 

We are now in a position to estimate the necessary 
values of the interaction parameters. Let us express the 
time-dependent Rabi frequency of the pulse as fl(t) = 
£lof(t/T), where fio is its peak value, T is the char- 
acteristic pulse width, and f{t/T) describes the pulse 
shape. The AE approximation demands |A| 3> Oo. The 
value of the phase is in general of order 0(1); this implies 
QlT ~ |A|. Hence we must have (fl T) 2 ~ |A|T > fi T, 
which in turn implies 



|A| » n Q » i/T, 



(14) 



i.e. the detuning A and the peak Rabi frequency Qq must 
be large compared to the Fourier bandwidth of the pulse 
1/T. 



C. Examples 
Gaussian pulse. For a Gaussian pulse, 

= VL e~ t2 ' T \ (15) 

where T and Qq are positive constants, the Schrodinger 
equation cannot be solved exactly and only some approx- 
imations are known [10]. The AE approximation (13) 
gives 



nlTy/n fi§0F(OgT 2 - 4v/2) 
4V2A 32A 3 T 



+O(n 6 /A 5 ). (16) 



Due to condition (14), the "roughness" term — 4\/2 in 
the second term, which derives from the fi(i)-term in 
Eq. (13), is small compared to the term f^T 2 . 

A meaningful phase gate requires that <j> ~ w, which 
implies that we must have QqT ~ 4V27r|A|. Then the 
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FIG. 2: (Color online) Phase shift vs the detuning A for 
a Gaussian pulse shape with a peak Rabi frequency flo = 
8/T. The exact phase, calculated numerically by solving the 
Schrodinger equation, is compared with the AE and AE2 ap- 
proximations, the adiabatic (AA) and superadiabatic (SA) 
phases. Top: the phase shift; bottom: the absolute error of 
the respective approximation. 



second term in Eq. (16) is ~ 7r 3 / 2 /(AT), which provides 
an estimate of the error in the phase. For an error < 
10~ 4 , we must have |A|T > 5 x 10 4 . This is indeed 
a very large value, particularly in many-particle systems 
where a variety of modes exists, and such a detuning may 
violate the condition of single-mode coupling. 

Figure 2 shows the phase shift after an interaction with 
an off-resonant Gaussian pulse. The AE phase (16) ap- 
proaches the exact phase only when the detuning A ex- 
ceeds the peak Rabi frequency fio The error of the first 
(steady-state) term in the AE phase barely drops to 1% 
in the shown range. The second term in the AE expan- 
sion (16) is seen to improve the accuracy as | A| increases. 
The other, more accurate phases shown in the same fig- 
ure are derived in the following section. 

Figure 3 shows the same phases versus the peak Rabi 
frequency flo. Similar conclusions can be drawn, as for 
Fig. 2: the AE approximation gives reasonable results 
only when the detuning A greatly exceeds the peak Rabi 
frequency f2o- The other two phases, to be discussed 
below, clearly provide much better fits to the exact phase 
for the entire ranges in Figs. 2-3. 

Hyperbolic secant pulse. The sech pulse, 



Q(t) = fl sech(t/T), 



(17) 



describes the pulse shape in the famous exactly-soluble 
Rosen-Zener model [111; we shall return to it in Sec. VI. 
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FIG. 3: (Color online) Phase shift vs the peak Rabi fre- 
quency Qo for a Gaussian pulse shape and a detuning A = 
10/T. The exact phase, calculated numerically by solving 
the Schrodinger equation, is compared with the AE and AE2 
approximations, the adiabatic (AA) and superadiabatic (SA) 
phases. Top: the phase shift; bottom: the absolute error of 
the respective approximation. 



The AE approximation (13) gives for it 
n 2 (tt 2 T 2 - 2 



iscch 



2A 



12A 3 T 



O{Q b /A b 



(18) 



This approximation is compared in Fig. 4 with the ex- 
act values. Similar conclusions as for the Gaussian pulse 
in Figs. 2 and 3 apply: the AE approximation provides 
a reasonable estimate for the phase shift only when the 
detuning A greatly exceeds the peak Rabi frequency flo- 
Keeping more terms in the expansion (18) improves the 
accuracy to some extent for large detunings, but this nei- 
ther extends the range of validity of this approximation 
nor reaches the accuracy of the adiabatic and superadi- 
abatic approximations, to which we turn our attention 
now. 



IV. ADIABATIC APPROXIMATION 

We shall now demonstrate that the adiabatic approxi- 
mation is a very accurate tool for calculation of the phase 
shift, with a vast domain of validity. 



A. Two-state adiabatic solution 

The adiabatic states <p+(t) and (fi-(t) are the 
eigenstates of the time-dependent Hamiltonian in the 



FIG. 4: (Color online) The same as Fig. 2 but for 
hyperbolic-secant pulse shape. 



Schrodinger representation (7), with eigenvalues 

A±(t) = i[A±A(t)]. (19) 



with X(t) = ^il(t) 2 + A 2 . The amplitudes in the adia- 
batic basis a(t) = [a+(t), a_ (t)) T are connected with the 
diabatic ones b(t) via the rotation matrix 



R(ft) 



cos ft sin ( 
— sin 8 cos ( 



(20) 



as b(t) = R(0(i))a(i), where ft = iarctan(fi/A). The 
Schrodinger equation in the adiabatic basis reads 



ihd t &(t) = n a (t)&(t), 



where 



H a (i) 



A_(t) -i6{t) 
i6{t) \+{t) 



(21) 



(22) 



If \6(t)\ < A+(t) - A_(t) = X(t), then the evolution is 
adiabatic and the propagator in the adiabatic basis reads 



U a = 



-iA_ 





,— 1A4 



(23) 



where A± = / ' X±(t) dt. We find readily that the prop- 
agator U in the original basis (2) is the phase gate (1), 
U = F, with the phase <fi = — A_, or explicitly, 



Wn{t) 2 + a 2 - A] dt. 



(24) 



This phase reduces to the AE approximation (8) for large 
detuning, |A| ^> f2 - However, Eq. (24) is valid also for 
|A| < Oq, provided the adiabatic approximation holds. 
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Figures 2-4 show that the adiabatic phase (24) pro- 
vides a considerable improvement of accuracy over the 
AE phase for all detunings |A| > 1/T. Unless a very 
high accuracy is required (error < 1CP 4 ) the adiabatic 
phase must suffice in applications. The key to the under- 
standing of the reason for its accuracy is hidden in the 
adiabatic condition. 



B. Adiabatic condition 

For adiabatic evolution, the nonadiabatic coupling 
\0(t)\ must be small compared to the splitting X(t), in or- 
der to suppress transitions between the adiabatic states. 
For a constant detuning, the adiabatic condition reads 



|f2(i)A| <c 2[fl(t) 2 + A 2 ] 3 / 2 . 



(25) 



For Gaussian [10] and sech [11] pulse shapes this condi- 
tion reduces to 



Gaussian: 
sech : 



|A| > A 



|A| > A = 



3V3T 
1 

3V6T 



y/HtoaT), (26a) 
(26b) 



Adiabatic evolution is achieved for a sufficiently large 
detuning. The sech pulse is obviously more adiabatic for 
it requires a lower detuning. The Gaussian pulse is less 
adiabatic; moreover, unlike the sech pulse it exhibits a 
logarithmic power broadening. 

The important message for the present context is that 
the adiabatic approximation requires a much lower value 
of the detuning than the AE approximation. Conditions 
(26) are only indicative: a more thourough analysis shows 
that the nonadiabatic deviation vanishes exponentially 
with the detuning [10-12]. Consequently, the necessary 
detuning increases logarithmically with the required ac- 
curacy. The implication is that the adiabatic phase gate 
can operate also at intermediate detunings, moreover, re- 
gardless of the value of ri , because the adiabatic condi- 
tion does not depend (or depends very weakly) on f^o- 
Thus the condition for the adiabatic phase gate is 



|A|T> 1; 



(27) 



then the effect of coherent population return [13-15] - 
the adiabatic return of the population to the initial state 
in the absence of a level crossing — ensures a negligibly 
small transition probability in the end of the interaction. 

We point out, however, that the perturbative estimate 
for the transient excitation is still 



n(t) 2 



2[Sl(t) 2 + A 2 ] 



(28) 



hence there may be a significant transient excitation un- 
less |A| ^ £lo- The adiabatic phase gate, therefore, can 
be used for Q > |A| only if the relaxation times of the 
two qubit states are large compared to the pulse duration 
- a condition that must be fulfullcd for any practical 
qubit for all types of operations. 



C. Superadiabatic phase 

If the evolution is not perfectly adiabatic, we can diago- 
nalize the adiabatic Hamiltonian (22) by the transforma- 
tion a(t) = K(x(t))s(t), where x(t) = \ arctan[0(t)/A(t)] 
and the vector s(t) contains the amplitudes in the su- 
peradiabatic basis; they satisfy the equation 



ihd t s(t) = U s (t)s(t), 



where 



H s (i) = h 



**(*) 



-»*(*) 
/*+(*) 



(29) 



(30) 



and (j,± = (A ± V^ 2 + A 2 +40 2 )/2. The condition for 
superadiabatic evolution is |x(t)| <C ^+{t) — fJ>-(t); if it 
holds then the propagator in the superadiabatic basis is 



-iM_ 







-iM + 



(31) 



where M± = f* ! fi±(t)dt. The propagator in the orig- 
inal basis (2) is the phase gate (1), U = F, with the 
superadiabatic phase (j) — —M-, or explicitly, 



' fl(t) 2 + A 2 + 



n(t) 2 A 2 

[n(t) 2 + A 2 ] 2 



dt. 



(32) 

The derivative term is the superadiabatic correction to 
the adiabatic phase (24). The condition this correction 
to be small is the same as the adiabatic condition (25). 

As evident from Figs. 2-4, the superadiabatic phase 
(32) is extremely accurate for all detunings except A — > 
0, with an error comfortably below 10~ 4 (the usual fault 
tolerance in quantum computing). The results clearly 
demonstrate that the superadiabatic phase outperforms 
the adiabatic phase, let alone the AE phase. 

Following the same diagonalization procedure, one can 
go to the next superadiabatic bases and achieve an even 
higher accuracy. However, as Figs. 2-4 suggest, this is 
unnecessary since the superadiabatic phase already easily 
satisfies the commonly accepted accuracy goal. 



D. Dissipation effects 

Dissipation is detrimental for quantum information 
processing and various proposals have been put forward 
to reduce its effects [1]. Because it is impossible to treat 
here all aspects of dissipation we restrict ourselves only 
to the simplest case of irreversible population loss. In 
femtosecond physics, where the present results can be 
particularly useful, this is a very reasonable assumption 
because population loss can occur through ionization in- 
duced by the driving laser pulse, whereas the other types 
of dissipation (dephasing, spontaneous emission, etc.) 
are irrelevant due to the ultrashort time scale. 
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In order to account for the population loss, we write 
the Hamiltonian in the form 



0(f) 
Q(t) 2A(t)-iT(t) 



(33) 



Next, we go to the adiabatic basis, formed of the eigen- 
states of the Hamiltonian without losses (r = 0) . In this 
basis the Hamiltonian reads 



H 



A_-iirsin 2 6> \irsm20- it 
\iTsin2e + ie A 



+ 



\iT cos 2 9 



(34) 



where for brevity the argument f is omitted. In the adi- 
abatic limit, we can neglect 9, as was done in Sec. IV A. 
Next, we recognize that the first (second) adiabatic state 
coincides at t — > ±oo with the first (second) diabatic 
state. We eliminate adiabatically the second adiabatic 
state and obtain the phase and the population of state 1, 



Pi = 



exp 



A_(f) 
tf 



T(t) 2 sin 2 29 (t) 
+ 16A(t) 

T(t) sin 2 9{t) dt 



dt, (35a) 
(35b) 



We note that the population is much more sensitive (ex- 
ponentially) to losses than the phase, which in the lowest 
order is quadratic in T. This conclusion is demonstrated 
in Fig. 5. The phase barely changes its value as the loss 
rate changes from zero to 10/T, even as the population 
of the initial state decreases considerably. We note that 
the displayed range of loss rates is much larger than what 
can be tolerated in quantum computing; nevertheless the 
phase itself and the approximation to it are very stable 
against such losses. 



EFFECTS OF ADDITIONAL STATES 



A. Ladder configuration 
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FIG. 5: (Color online) The initial-state population and phase 
vs the loss rate Y for a Gaussian pulse with Qo = 8/T and 
A = 20/T. The solid lines represent the exact values and the 
dashed lines are for the approximations (35). 



As expected, in the limits fl 23 — > or |A 3 | 
expression reduces to Eq. (8). When 4|A 2 A 3 | 
we obtain 



-> oo, this 

» o 23 (t) 2 



4A 2 



dt 



tf n 12 {tfn 23 {tf 

16A 2 A 3 



dt 



(38) 



We conclude that for 0x2 ~ ^23 and A 2 ~ A3, the cor- 
rection from the presence of an additional state is of the 
same order as the second term in the AE expansion (13). 

Adiabatic approximation. For simplicity and with- 
out loss of generality, we will assume A 3 > A 2 > 0. In 
order to find the eigenvalues of the Hamiltonian (36), we 
have to solve the cubic characteristic equation 



Adiabatic elimination. We shall now find how the 
presence of an additional state affects the phase shift. For 
this purpose let us consider a three-state system in the 
ladder configuration, wherein nonzero dipole moments 
only link state 1 with state 2 and state 2 with state 3 
(Fig. 1, middle), as described by the Hamiltonian 



Ox 2 (t) 

n 12 (t) 2A 2 n 23 (t) 
n 23 (t) 2 a 3 



(36) 



For the first-order (steady-state) AE approximation, we 
set c 2 (t) = and 03(f) = 0, and find from the Schrodinger 
equation the accumulated phase in state 1 to be 



A 3 0i 2 (t) 2 



4A 2 A 3 -0 23 (f) 2 



■dt. 



(37) 



as 



be + c = Q, 



(39) 



where a 

-)2 



-A 2 - A 3 , b = A 2 A 3 - (n\ 2 + 2 3 )/4, c - 
A3^i2/4- The three roots of this equation are the 
quasienergies of the three-state system [9] , 



£1 
£3 



where 



p = \J a 2 — 3b, 
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FIG. 6: Phase shift of the lowest state in a three-state ladder 
with detunings A3 = 2A2 = 20/T vs the peak Rabi fre- 
quency Qo- The Rabi frequencies for the two transitions are 
equal and have Gaussian shapes. The phases from the adia- 
batic approximation for two (AA2) and three (AA3) states, 
Eqs. (24) and (42), and the AE approximation for three states 
(38) are compared to the exact values. The AA3 curve is 
nearly indiscernible from the exact one. 



FIG. 7: Phase shift of the lowest state in a three-state V- 
system with detunings A3 = 2A2 = 20/T vs the peak Rabi 
frequency fio- The Rabi frequencies for the two transitions 
are equal and have Gaussian shapes. The phases from the adi- 
abatic approximation for two (AA2) and three (AA3) states, 
and the AE approximation for three states (45) are compared 
to the exact values. The AA3 curve is indiscernible from the 
exact one. 



B. V configuration 



The adiabatic phase is just an integral over £\{t) [24], 



L (t)dt. 



(42) 



In order to exhibit the effect of the third state on the 
phase gate, we derive the asymptotics of the quasienergy 
£1 for large A3 and substitute it in Eq. (42); we find 



it £ - {t)dt ~J t 



e-(t)n 23 (t) 2 



-.At, 

(43) 



where e_ 



(A 2 



ny/2. 



Figure 6 compares the AE and adiabatic approxima- 
tions to the exact phase shift in a ladder system. The 
three-state adiabatic phase (42) is very accurate through- 
out, as already anticipated, whereas the AE phase (38) 
departs from the exact one when the Rabi frequency be- 
comes comparable and exceeds the detunings. We note 
that the error of the AE phase is of the same order as the 
error from the neglect of the additional state, Eq. (24). 



If the three-state system is in a V configuration (Fig. 
1, right), the Hamiltonian reads 



n 12 (t) n 13 (t) 

fiia(t) 2A 2 

n 13 (t) 2A 3 



(44) 



Then the AE approximation, applied in a similar fashion 
as for the ladder system above, gives in the first order 
the expression 



4A 2 (i) 



dt- 



4A 3 (i) 



dt. 



(45) 



The contributions from each arm in the V system are 
independent. Higher-order terms mix the contributions 
from the two arms. 

The adiabatic phase is calculated in a similar manner 
as for the ladder system, in the form of an integral over 
the respective eigenenergy. 

In Fig. 7 we compare various expressions for the dy- 
namic phase in a three-state V system. The figure reveals 
that, once again, the three-state adiabatic approximation 
(AA3) provides a very accurate estimate. On the con- 
trary, the error of the AE approximation is comparable 
to the effect of the additional third state. 

To conclude this section, we point out that a third 
state in the considered ladder or V systems distorts the 
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symmetry of the phase gate (1) because the third state is 
coupled differently to the two qubit states: directly to one 
of them and indirectly (via a two-photon transition) to 
the other state. For each qubit state, the additional state 
will make the linkage look either as a ladder (if connected 
to the other qubit state) or V (if connected to the same 
qubit state). Consequently, the phase shifts for the two 
qubit states are different, as evident when one compares 
the expressions for the ladder and V systems above. We 
note that a Lambda-system, with the population initially 
in one of the lower states is equivalent in the present 
context to a ladder system (population initially in the 
end of the chain) . The phase shifts in such a system have 
been studied recently [8], in a slightly different context. 

In the next section we compare these approximations 
with the values of several exactly-soluble two-state and 
three-state models. 



VI. EXACT SOLUTIONS 

We shall present the solutions of three exactly soluble 
models. The first one is the famous Rosen-Zencr (RZ) 
model, which assumes a sech pulse shape(17) and a con- 
stant detuning, and the others extend the RZ model to 
systems with three states. The latter are solved by using 
the Majorana decomposition [16, 17] for a ladder system 
and the Morris-Shore transformation [18] for a V system. 



A. Two states: Rosen-Zener model 

The exact solution for the phase in the RZ model (17) 
is [11, 19] 



arg 



r 



, (46) 



where a — SIqT, S = AT, and T(z) stands for the gamma 
function [20]. Using the Stirling asymptotics of T(z) [20], 
we obtain 



<f>. 



2A 



^IT 2 2) 
12A 3 T 



+ ... 



|A|»1/T,Q ). (47) 



This is exactly the result from the AE expansion (13). 

The exact expression for the phase (46) allows us to 
perform a theoretically exact phase gate operation with 
a variable phase <j) by selecting a suitable detuning A. By 
using standard properties of the gamma functions, one 
can show that the transition probability vanishes exactly 
for a pulse area A — ttQqT = 2mr, where n is an integer. 
For these values (a = 2n), Eq. (46) reduces to [19] 



nil + 2 arg 



1 - iAT) 



k=l 



(48) 



For any desired phase (f> the corresponding detuning is 
found by solving the latter equation for A. For example, 



a phase shift <p — it can be obtained by any odd value 
of n = 1,3,5,... and A = (a property that is well 
known for any pulse shape). The same phase shift can 
also be produced by a pulse area A = 4ir (n = 2) and 
AT = \/3. As further examples, phase shifts of 4> = tt/2, 
7r/3, 7r/4, and tt/6 can be obtained for n = 1 by choosing, 
respectively, AT = 1, V3, 1 + V2, and 2 + V^. 



B. Three-state ladder 

We consider a three-state system in a ladder config- 
uration (Fig. 1, middle), described by the Hamiltonian 
(36), assuming that 

n 12 (t) = n 23 (t) = n sech(t), a 3 = 2A 2 = 2A. (49) 

This model has a simple exact solution [21]: the am- 
plitude of state 1 is just the square of the amplitude 
for a two-state problem with a Rabi frequency fl(t) = 
Qi2(t)/V2- Because this is the two-state RZ model, we 
find the phase of state 1 to be 



2 arg 



(50) 

with a = O T and 5 = AT. The exact expression (50) 
allows one to design an exact phase gate even in the pres- 
ence of an additional state. For example, a gate phase 
= 7r can be obtained for f2 T = 2ny/2 with any odd 
value of n = 1, 3, 5, . . . and A — 0; the same phase can 
be produced by a pulse area A — 4ir\/2 (n = 2) and 
AT = \/3. A phase shift of 4> = tt/2 can be obtained for 
A = 2irV2 and AT = 1. 

The asymptotics of this phase reads 



2A 



%{%T 2 4) 
24A 3 T 



+ . 



|A|»l/T,O )- (51) 



The comparison with Eq. (47) shows that the effect of 
the third state emerges only in the second term in the 
expansion over A, a feature that appeared earlier in the 
AE expansion (38) and in the adiabatic expansion (43). 



C. Three-state V-system 

If the system is in a V configuration (Fig. 1, right), 
and if 



fli2 = sech(t/T), 
^13 = >*i3^o sech(i/T), 
A 2 = A 3 = A, 



(52a) 
(52b) 
(52c) 



with and >fi 3 arbitrary constants, then we can use 
a simple change of basis, known as the Morris-Shore 
transformation [18]. The latter transforms the V-system 
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into an uncoupled dark state and a two-state system, 
with a detuning A and a Rabi frequency x£l{t), where 
k = \J x\ 2 + xf 3 . In this manner it can readily be shown 
that the argument of the amplitude of state 1 is [22] 



2 arg 



xa 



The asymptotics reads 



(53) 



2A 



x 2 n 2 {x 2 n 2 -2) 

12A 1 



+ ... (|A| > 1/T,fi ). 

The leading term is a sum of two independent phase shifts 
induced by each of the two arms of the V-system, in 
agreement with Eq. (45), while the ^ 4 -part in the second 
term represents a combined contribution of the two arms. 



VII. CONCLUSIONS 

We have presented a detailed analysis of the accuracy 
of the adiabatic elimination approach, which allows one 
to eliminate weakly coupled far-off-resonant states and 
reduce the interaction dynamics of a quantum system to 
a smaller effective one. We have put a special emphasis 
on the acquired phase shifts in the probability ampli- 
tudes after the interaction with an off-resonant pulsed 
field. The results have direct implications in the con- 
struction of variable dynamic phase gates, which are of 
major importance in quantum information processing. 

We have shown that in the traditional implementa- 
tion with a far off-resonance pulsed field, the formula for 
the gate phase derived by the adiabatic-elimination ap- 
proximation, has to be used with great care because it 
is just the first term of a series expansion in the inverse 
detuning 1/A; higher terms, unless negligible in value 
with respect to unity, may render the formula irrelevant. 
This formula requires a very large detuning in order to be 
sufficiently accurate, because including corrections from 
higher terms is barely useful due to their complexity and 
slow convergence. However, quantum information pro- 



cessing involves operations with entangled many-qubit 
systems, possessing a variety of frequency modes; there a 
very large detuning may violate the assumption of single- 
mode interactions. 

We have proposed to use a much more accurate formula 
for the phase shift, derived within the adiabatic-following 
approximation. The advantage of this adiabatic phase is 
that the condition for adiabatic evolution is much more 
relaxed than the condition for adiabatic elimination. The 
adiabatic phase contains the AE approximation as a lim- 
iting case for |A| — ► oo; however, it also applies to mod- 
erate detunings. The adiabatic phase applies also to the 
case when the peak Rabi frequency f2o exceeds the de- 
tuning A; in this case the phase gate operates due to 
the effect of adiabatic population return. A superadia- 
batic correction is demonstrated to further improve the 
accuracy, to errors comfortably below the fault-tolerance 
limit of 10~ 4 in quantum computing. We have also de- 
rived the corrections to the gate phase from additional 
states coupled to the two qubit states. 

In addition, we have derived several exact expressions 
for the gate phase in several exactly soluble analytical 
two-state and three-state models, assuming a hyperbolic- 
secant pulse shape. The exact analytical formulae allow 
us to design highly accurate phase gates, however, at the 
expense ot the requirement for a special pulse shape. We 
have also used the exact expressions to test the accuracy 
of the derived AE and adiabatic approximations. 

The results in this paper have potential applications 
not only in the calculation of dynamical Stark phase 
shifts in simple phase gates but also in complicated multi- 
state linkage patterns, which can be factorized to simpler 
systems by utilising the intrinsic symmetries, e.g. by the 
Morris-Shore transformation [18, 23]. 
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